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Active/Adaptive Flexible Motion Controls with Aeroservoelastic System Uncertainties 

Finite Element Model Tuning of X-56A Aircraft using Parallelized Big-Bang Big-Crunch Algorithm 
Unsteady Aerodynamic Model Tuning Based on In-direct Method 
Design of an Aeroservoelastically Tailored Wings and Aircraft 
Aeroelastically Tailored Wing Designs 
Aeroservoelastically Tailored Wing Designs 
Reduced Order Modeling 

Equivalent Beam Modeling for X-56A Flight Simulations 
Development of CFD based Flutter Analysis Technique 
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Active/Adaptive Flexible Motion Controls 
with Aeroservoelastic System Uncertainties 


Research Goals/Obiectives: 

□ Aeroservoelastic model validation is an essential 
procedure for the safety of flight. 

□ Uncertainties still exist in aeroservoelastic system 
even with the test validated aeroservoelastic model 
due to 

*t* * time-varying uncertain flight conditions, 

*t* transient and nonlinear unsteady aerodynamics and 
aeroelastic dynamic environments. 

□ For a flexible motion control problems, we need a 
control law that adapts itself to such changing 

conditions. 

Approach: 

□ Digital adaptive controller for a flexible motion control 
♦> On-line Parameter Estimation & Health Monitoring 
*t* On-line Control Law Design 

> Baseline Control Law Design will be based on the 
test validated aeroservoelastic model 

> On-line modification of control law will be based 
on the estimated "delta system model". 

□ Model Validation and Tuning 

*t* Minimizing uncertainties in an aeroservoelastic 
model 

Structural Dynamic Model Tuning 
Actuator Model Tuning 

> Unsteady Aerodynamic Model Tuning 

Team: Chan-ei Pak. Martv Brenner. Roger Truax. & Alex Chin 
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Finite Element Model Tuning of X-56A Aircraft 
using Parallelized Big-Bang Big-Crunch Algorithm 



Objectives 


□ The primary objective of this study is to reduce uncertainties in the structural dynamic finite element model of an 
aircraft to increase the safety of flight. 

□ This model tuning technique is applied to improve the flutter prediction of the X-56A aircraft. 

□ This work is supported by the Aeronautics Research Mission Directorate (ARMD) Fixed Wing and High Speed 
projects under Fundamental Aeronautics (FA) program. 



□ Collaboration with AFRL & LMSW 
<♦ Two Center Bodies 
*X* One Rigid Wing 
<♦ Three Flexible Wings 
Ground Control Station 
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X-56A Aircraft: Exploded View of X-56A 


Bay 

Volume 

(ft^) 

% 

Occupied 

Comment 

A 

1.90 

10% 

Nose Wheel System Bay 

Aa 

20% 

Power System Bay 

B 

1.10 

0% 


Bb 

0% 


C 

2.00 

70% 

Electronics Bay 

Cc 

90% 

Ballistic Recovery System 

D 

0.46 

30% 

Fuel System 

Dd 

30% 

Fuel System 

E 

0.94 

100% 

Fuel Wet Bays 

Ee 

100% 


6.9 ft^ Useable volume (54% of total) 



Wheel pants for 
increased speed 


Wing Attach Fitting 



Forward and aft ballast bays for stability tuning 


12 lb Water 
ballast each 

611b Water 
ballast 

Wing Cameras 
Nose Camera 
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Flutter Analysis Procedure @ NASA Dryden 


Everyone believes the test data except for the experimentalist, and no one believes the 
finite element model except for the analyst. 




Some of the discrepancies come from analytical Finite Element modeling 
uncertainties, noise in the test results, and/or inadequate sensor and actuator 
locations. Not the same orientation for each sensor. 


Structural Dynamic 
Finite Element Model 


Weight, C.G., Moment of 
inertia, & GVT data 


□ Flutter Analysis 






Uncertainties in the structural dynamic model are minimized through the use of 
"model tuning technique" 


□ 


Based on analytical modes 
Validate Structural Dynamic Finite Element Model using Test Data and Update if needed 




Use MDAO (Multidisciplinary Design, Analysis, and Optimization) tool with Model 
Tuning Capability or Standalone Model Tuning Code 

> Model tuning is based on optimization. 

^ Design Variables 

• Structural sizing information: Thickness, cross sectional area, 
area moment of inertia, etc. 

• Point properties: lumped mass, spring constant, etc. 

• Material properties: density. Young's modulus, etc. 

Constraints 


Structural Dynamic 
Model Tuning 


Validated Structural 
Dynamic Model 


Create Unsteady 
Aerodynamic Model 
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□ Approach 


Structural Dynamic Model Tuning 
using Object Oriented Optimization Tool 




□ 


Minimize "objective functions" using object oriented optimization (03) tool 
which leverages existing tools and practices, and allows 
the easy integration and adoption of new state-of-the-art 
software. 

Optimization Problem Statements 


Minimize J = 


Optimizer 


i ;/ 


G(x) 


Optimization 
Tool 


> /: Objective function 

> w,: Weighting factor 
for the performance 
index i 

> Jf. Performance index i selected 
for objective function 

> Jj^. Performance index ^selected for constraint functions 

> Sj^. Small tolerance value for performance index k 
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Measure Weight, 
C.G., Moment of 
inertia. 

Frequencies, & 
Mode Shapes 
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Parallelized Big Bang Big Crunch Algorithm 




□ A global optimizer 

First step: Big Bang step 

> Selection of the N (number of population) random design variable vectors X, 
(i=2,3,...,N) using uniform random number generator such that 

XLi<Xi<XU, 

> Current design configuration is saved in the design variable vector X^. 

> Create M ( N = M x integer ) number of design variable vectors simultaneously. 

> M number of objective functions will be computed simultaneously. 

Needs M number of NASTRAN licenses 

Second step: Big Crunch step 

> Shrink design variable vectors to a single representative design point via a center 




of gravity (CG) 


^CG - 


i=\ 

1 




2 j 

i=\ 'Ji 

Third step: Big Bang step 

> Compute new candidate design variable vectors around the CG location using the 
standard normal random number generator 

m(XU,- - XL. ) 


x^=/?Xcg+(1-A)Xgo + - 


NBB 


where, r is the standard normal random number; a is the parameter limiting 
the size of the design space; NBB is the number of current big bang iteration; 
and 3 is the parameter controlling the influence of the global optimum 
solution Xqo- 

Parameters a and 3 for the best performance was a=l and 3=0.2 for the truss 
design problems and a=l and 3=0.7 for the parameter estimation problems. 
Go to the second step until converge 
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Pre-GVT Modal Analysis of X-56A 



structural Dynamic Finite Element Model 


□ Based on MSC/NASTRAN code 
Assembled configuration 
<♦ 8249 nodes 

Use 40 modes for the flutter analysis 
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Modal Analyses of Each Structural Component 



Main Landing Gear 


(a) Right Wing: 3325 nodes 


Nose Landing Gear 
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Mode Shapes 



Cantilevered Wing 


Mode 9: 9.723 Hz 


Mode 7: 3.499 Hz 


Mode 1: 2.298 Hz 


Modes: 5.557 Hz 


Mode 2: 9.669 Hz 


Mode 10: 10.94 Hz 
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Modes: 18.58 Hz 


Mode 11: 11.94 Hz 


Mode 27: 38.26 Hz 


Mode 15: 15.16 Hz 


Mode 31: 48.51 Hz 
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Mode Shapes (continued) 


Mode 6: 31.70 Hz 


Mode 30: 43.C 


Mode 7: 38.81 Hz 


Mode 8: 48.84 Hz 


Mode 


15.43 


Hz 


Mode 4: 14.82 Hz 


lode 3: 11.48 Hz 


Mode Shapes: Boom & Landing Gears 



Mode 1: 13.51 Hz 


Mode 2: 15.36 Hz 


Mode 3: 15.43 Hz 
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Pre-GVT Flutter Analysis of X-56A 



Unsteady Aerodynamic Model 



Flow 


□ Based on ZAERO code 

<♦ 416 elements 

Select 16 reduced frequencies between 0 & 1 
Mach = .130, .195, and .284 
*X* Linear Theory 
*t* Use Matched Flutter Analysis 


Splining Points 
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Frequency (Hz) Damping 


V-g and V-f Curves at Mach = 0.130 Before Model Tuning 



Dampin 


Modal Participation Factors 


Mode 


Freq. 

(Hz) 


Flutter 


Flutter 


Flutter 


0.000 


30.3 


73.3 


25.9 


87.1 


40.9 


40.9 


3.499 


15.8 


89.1 


94.7 


5.557 


26.1 


67.0 


9.723 


43.0 


43.0 


61.2 


61.2 


10.94 


26.1 


93.1 


11.94 


97.7 


97.2 


15.16 


94.4 


96.1 


38.26 


92.1 


99.2 


43.00 


95.4 


48.51 


95.6 


98.2 


19.01 


96.7 


Flutter Mode 

Speed 

Frequency 

Altitude 

1 

101.0 Keas 

2.597 Hz 

-9204.3 ft 

2 

110.1 Keas 

7.134 Hz 

-14330. ft 

3 

121.0 Keas 

4.971 Hz 

-20169. ft 
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First Flutter Mode Shapes of X-56A Aircraft 
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Second Flutter Mode Shapes of X-56A Aircraft 
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Third Flutter Mode Shapes of X-56A Aircraft 
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Altitude (ft) 


Flutter Boundaries 
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Updating the Finite Element Model of the X-56A 
Using Ground Vibration Test Data 



structural Dynamic Model Tuning using GVT Data 


□ NASA Standard: NASA-STD-5002 Section 4.2.6.d 

♦> NASA Technical Load Analysis of Spacecraft and Payloads 

*t* Agreement between test and analysis natural frequencies shall, as a goal, be within ^ for the significant 
modes. 

*t* Accurate mass representation of the test article shall be demonstrated with orthogonality checks using the 
analytical mass matrix M and the test mode shapes Oq. The orthogonality matrix is computed as Oq’'’M Og. As 
a goal, the off-diagonal terms of the orthogonality matrix should be less than ^ for significant modes based 
on the diagonal terms normalized to 1.0. 

□ Military Standard: MIL-STD-1540C Section 6.2.10 

♦> Test Requirements for Launch, Upper-Stage, & Space Vehicles 

*t* Analytical model frequencies are to be within 3% of test frequencies. 

*t* Using a cross-orthogonality matrix formed from the analytical mass matrix and the analytical and test modes, 
corresponding modes are to exhibit at least 95% correlation and dissimilar modes are to be orthogonal to 
within 10% . 

□ AFFTC-TIH-90-001 (Structures Flight Test Flandbook) 

♦> If measured mode shapes are going to be associated with a finite element model of the structure, it will 
probably need to be adjusted to match the lumped mass modeling of the analysis. 

*t* Based on the measured mode shape matrix Og and the analytical mass matrix M, the following operation is 
performed: 3>g’'’M Og 

♦> The results is near diagonalization of the resulting matrix with values close to 1 on the diagonal and values 
close to zero in the off-diagonal terms. Experimental reality dictates that the data will not produce exact unity 
or null values, so 10 percent of these targets are accepted as good orthogonality and the data can be 
confidently correlated with the finite element model. 
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Bottom-Up Approach for Structural Dynamic Model Tuning 


Cantilevered or Free-Free Right Wing with Winelet 



Constrained Center-Bodv 

Front View 
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Use lyipdel Tuning 


X-56A Aircraft with Soft Suspension System 

Hoists 




Assembled X-56A 
Testbed Aircraft 


Use Model Tuning 


Use Mirror Image 
Use Model Tuning' 



Center Body Model 











- 
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Structural Dynamic Model Tuning 


□ Nose boom 

Beam elements 

Match first two modes, 15.36 Hz and 15.43 Hz; improve CG; and mass matrix 

□ Main landing gear 

*t* Beam elements 

*t* Match first two modes, 13.51 Hz and 16.75 Hz ; improve CG; and mass matrix 

□ Nose landing gear 

*t* Beam elements 

*t* Match first two modes, 33.03 Hz and 33.08 Hz ; improve CG; and mass matrix 

□ Center body 

*t* Improve CG and mass matrix 

□ Wing with winglet 

*t* Match first eight modes 

^ Use different weighting factors 
^ Modes 1 & 2 : 1.0 

^ Modes 4, 7, & 8 : 0.5 
Modes 3, 5, & 6 : 0.01 
Most time consuming part 

□ Assembled configuration 

Use superelement 

*t* Match 6DOF springs between wings and center body 
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Flutter Analysis Using Validated Structural 
Dynamic Finite Element Model 




Unsteady Aerodynamic Model Tuning 
Based on In-direct Method 


□ 






Unsteady Aerodynamic Model Tuning based on Direct Method 
The NASA Dryden has developed an Object-Oriented Optimization (0^) tool. 

The 0^ tool leverages existing tools and practices, and allows the easy integration and adoption of new state-of- 
the-art software. 

Local gradient based optimizer as well as 
global optimizers are available. Hybrid methods 
are also available. 




Optimizers: 

DOT (local). Genetic Algorithm (GA), & 
Big Bang-Big Crunch (BBBC) algorithm 




Hybrid optimizers: 
GA(CDV)+DOT(CDV), 
GA(CDV)+DOT(CDV)+GA(DDV)) 
BBBC(CDV)+DOT(CDV), & 
BBBC(CDV)+DOT(CDV)+BBBC(DDV 


□ In-direct Method 

Change AlC through the change of aerodynamic 
panel geometry 




Frequency 

Difference 



Frequencies 
Measured 
from Flight 
Test 
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Design of an Aeroservoelastically Tailored Wings and Aircraft 


Research Goals/Obiectives: 

□ Use aeroelastic tailoring theory and active flexible 
motion control technique to satisfy the overall strain and 
aeroelastic and aeroservoelastic instability requirements 
within given flight envelopes 

□ Use curvilinear sparib concept as well as composite ply 



Approach: 

□ Simultaneously update structural as well as control 


design variables during early design phase 
*X* Perform topology optimization with curvilinear sparibs 
*1* Use aeroelastic tailoring up to Vd line 
*t* Use aeroservoelastic tailoring between Vd and 1.15 Vd 

□ During MDAO code development use X-56A as a sample 
case 

□ Design ARIO Wing using object-oriented MDAO tool 

*1* Design scaled ARIO wing using structural model tuning 
tool 

*1* Candidate aircraft: B-777 or B-747 (SOFIA) 

□ Design ARM Wing using Object-Oriented MDAO tool 

*t* Design scaled ARM wing using structural model tuning 
tool 

<♦ Candidate aircraft: Truss Braced Wing, etc. 

□ Design Supersonic Biplane 

Structural Dynamics Group 




Team: 


High Speed 


Aeroelastic stability envelope 

Flutter Boundary 




Vo 1-15V 

■ 1.0 
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Use 

Aeroelastic 
Tailoring 
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Aeroservoelastic 

Tailoring 


□ Chan-gi Pak, Wesley Li, Alex Chin, Roger Truax, & Paul Yoo 
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Equivalent Beam Modeling for X-56A Flight Simulations 


Research Goals/ Objectives: 

□ Create accurate and affordable model for 

aeroservoelastic simulations 

Approach: 

□ Create a wing equivalent beam or plate model using 
structural dynamic model tuning tool 

♦> Use equivalent beam model for high aspect ratio 
wing 

Use equivalent plate model for low aspect ratio wing 

Applications: 

□ Wing equivalent beam modeling 
<♦ High aspect ratio wings 

> X-56A 

> ARlOWing 

> ARM Wing 

□ Wing equivalent plate modeling 
♦> Low aspect ratio wings 


Reduced Order Structural Dynamic Modeling 
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Team: 

□ Chan-gi Pak & Peter Suh (Ph.D. Student @ Georgia Tech) 
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Development of CFD based Flutter Analysis Technique 


Research Goals/ Objectives: 

□ A flutter analysis technique in the transonic flight regime. 
The technique uses an iterative approach to determine the 
critical dynamic pressure, i.e. flutter boundary, for a given 


i/^ach number. 

kAA I 

h A A'"' 

I ' 

Unlike other CFD- 

rV V * 

based flutter analys 

fVV' * 

is methods, each 


iteration solves for the critical dynamic pressure and uses 
this value in subsequent iterations until the value 
converges. This process reduces the iterations required to 
determine the critical dynamic pressure. To improve the 
accuracy of the analysis, the technique employs a known 
structural model, leaving only the aerodynamic model as 
the unknown. 



Aporoach: 

□ The known structural model is represented as a FEM. 

□ The unsteady CFD analysis is performed. The output time 
history of the surface pressure is converted to a nodal 
aerodynamic force vector. The forces are then normalized 
by the given dynamic pressure. 

□ A multi-input multi-output parameter estimation 
software. Eigenvalue Realization Algorithm, estimates the 
aerodynamic model. 

□ The critical dynamic pressure is then calculated using the 
known structural model and the estimated aerodynamic 
model. 

□ This output is used as the dynamic pressure in subsequent 
iterations until the critical dynamic pressure is 
determined. 

Team: 

□ Chan-gi Pak & Paul Yoo (Ph.D. Student @ DSC) 
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Local Mach number 
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Questions? 


IG Trim Analysis 
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Frequency (Hz) Damping 


0.1 


V-g and V-f Curves at Mach = 0.195 Before Model Tuning 



1 % 
Damping 


0 20 40 60 80 


100 120 140 160 180 200 


Specd(Keas) 



Specd(Keas) 


Mode 

Frequency 

1st 

2nd 

3rd 

1-6 

0.000 Hz 

2 

28.5 

2 

27.1 

1 

31.6 

7 

3.499 

3 

13.4 

3 

6.8 


0.0 

8 

5.557 


0.0 


0.0 

2 

31.3 

9 

9.723 

1 

49.4 

1 

61.7 


0.0 

10 

10.94 


0.0 


0.0 

3 

31.3 

11 

11.94 

8 

0.8 

6 

0.9 


0.0 

15 

15.16 

5 

1.7 

5 

0.9 


0.0 

20 

19.01 

7 

0.9 


0.0 


0.0 

27 

38.26 

4 

2.4 

7 

0.8 


0.0 

30 

43.00 


0.0 


0.0 

4 

2.2 

31 

48.51 

6 

1.0 

4 

1.0 


0.0 


Flutter Mode 

Speed 

Frequency 

Altitude 

1 

94.87 Keas 

2.553 Hz 

16045. ft 

2 

110.8 Keas 

6.814 Hz 

8153.8 ft 

3 

111.9 Keas 

4.807 Hz 

7636.9 ft 
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Frequency (Hz) Damping 


V-g and V-f Curves at Mach = 0.284 Before Model Tuning 



Mode 

Frequency 

1st 

2nd 

3rd 

1-6 

0.000 Hz 

2 

27.8 

2 

28.3 

3 

24.7 

7 

3.499 

3 

11.4 

3 

6.1 


0.0 

8 

5.557 


0.0 


0.0 

1 

34.9 

9 

9.723 

1 

53.2 

1 

62.0 


0.0 

10 

10.94 


0.0 


0.0 

2 

34.9 

11 

11.94 

8 

0.7 

5 

0.7 


0.0 

15 

15.16 

5 

1.4 

7 

0.7 


0.0 

20 

19.01 

7 

0.8 


0.0 


0.0 

27 

38.26 

4 

2.1 

6 

0.7 


0.0 

30 

43.00 


0.0 


0.0 

4 

2.3 

31 

48.51 

6 

0.9 

4 

0.9 


0.0 


Flutter Mode 

Speed 

Frequency 

Altitude 

1 

91.67 Keas 

2.516 Hz 

34749. ft 

2 

111.7 Keas 

6.519 Hz 

26092. Ft 

3 

108.7 Keas 

4.725 Hz 

27354. ft 
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